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ABSTRACT 

This paper describes the application of linear 
stochastic optimal control theory to the design of 
the control system for the air intake (inlet) of a 
supersonic air-breathing propulsion system. The 
controls must maintain a stable inlet shock posi- 
tion in the presence of random airflow disturbances 
and prevent Inlet unstart, Two different linear 
time Invariant control systems are developed. One 
is designed to minimize a nonquadratic index, the 
expected frequency of inlet unstart, and the other 
is designed to minimize the mean square value of 
inlet shock motion. The quadratic equivalence 
principle is used to obtain the best linear con- 
troller that minimizes the nonquadratic performance 
index. The two systems are compared on the basis 
of unstart prevention, control effort requirements, 
and sensitivity to parameter variations. It is 
concluded that while controls designed to minimize 
unstarts are desirable in that the index minimized 
is a physically meaningful quantity, computation 
time required is longer than for the minimum mean 
square shock position approach. In addition, the 
simpler minimum mean square shock position solution 
produced expected unstart frequency figures which 
were not significantly worse than those of the non- 
quadratic solution, 

I INTRODUCTION 

The purpose of this paper is to study two ap- 
proaches to supersonic inlet control system design. 
Both approaches are similar in that they use opti- 
mal control theory, however each uses a different 
performance criterion. The overall goal is to 
provide a more rational and unified approach to in- 
let control design. 

A supersonic inlet is that portion of a supersonic 
propulsion system which decelerates air from 
supersonic velocity (ahead of the aircraft) to 
subsonic velocity (at the entrance to the com- 
pressor) . This deceleration is needed because 
present compressors require air at subsonic ve- 
locity. An inlet is a critical part of a super- 
sonic propulsion system. This is not the case for 
an inlet on a subsonic aircraft. This is because 
the dynamic head is large at high supersonic Mach 


numbers and may comprise a large percentage of the 
overall engine compression. For subsonic propul- 
sion systems, however, almost all the compression 
is done by the engine’s compressor. To aid the 
supersonic inlet in operating at peak efficiency 
in the face of varying flight conditions, certain 
variable geometry features, and associated con- 
trols, are required. 

A typical axisymmetric supersonic inlet is shown 
in Fig. 1 in a normal operating configuration. Air 
at supersonic velocity enters the Inlet past a weak 
oblique shock wave. It is compressed supersonical- 
ly past a minimum area point, or throat, up to the 
terminal normal shock. Thereafter, the flow is 
subsonic up to the compressor face station. 

A stable operating condition for the inlet is one 
in which the throat Mach number is greater than 
one and the normal shock is downstream of the 
throat. This is the so-called started condition. 

An upstream or downstream flow disturbance may, 
however, cause the throat Mach number to drop to 
one, or it may cause the normal shock to move 
ahead of the throat. When either of these occur, 
the inlet unstarts and enters an undesirable, un- 
stable operating region (called unstart). 

During an unstart a shock wave sweeps out of the 
throat and a strong shock wave forms ahead of the 
inlet. The resuLt is a large increase in drag and 
a large decrease in the pressure recovered at the 
compressor face. In addition there may exist an 
oscillatory flow pattern within the inlet . Air- 
craft performance under these conditions is poor 
and often unacceptable. Thus, unstart is very un- 
desirable. 

In order to prevent unstart and to maintain the de- 
sired operation of the inlet, two basic control 
modes are required. The first, shown in Fig. 1, is 
a translating centerbody, A forward translation 
of the centerbody causes the throat area to in- 
crease. By varying the centerbody position, the 
throat area can be varied so that the throat Mach 
number can be kept above one, even if the free 
stream Mach number decreases. The second manipu- 
lated variable is bypass door opening. Bypass 
doors are controlled to dump excess air in case a 
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disturbance causes a pressure rise at the compres- 
sor face. By spilling excess air, the pressure 
disturbance is prevented from pushing the normal 
shock forward ahead of the throat. 

Primary inlet disturbances can occur either up- 
stream or downstream of the inlet. Upstream dis- 
turbances may result from such things as shock 
waves from passing aircraft or atmospheric turbu- 
lence. Turbulence may consist of pressure, tem- 
perature or velocity (gusts) changes. Downstream 
disturbances may be due to changes in engine air 
flow demand which results from a pilot induced 
throttle change. Also aerodynamic noise may be 
generated by the compressor, or combustion noise 
may get fed back to the compressor face in duct 
burning fan type engines. The overall control de- 
sign goals are; (l) to keep the throat Mach num- 
ber and shock position as close as possible to 
their unstart limits. This maximizes efficiency. 
(.2) to keep the deviations of throat Mach number 
and shock position small so as to prevent unstart. 
(3) to design a control which produces desired re- 
sults while using a reasonable amount of control 
power. 

II PRESENT INLET CONTROL DESIGN TECHNIQUES 


The following sections of the paper are outlined 
as follows: Section III discusses the particular, 

somewhat restricted, optimal control problem con- 
sidered. The next section is an outline of linear 
quadratic optimal control and estimation theory 
and the application of this theory to nonquadratic 
performance indexes. Section V describes the 
computer algorithms used followed by section VI in 
which state variable inlet models and disturbance 
descriptions are presented. The results of a con- 
trol study comparing and contrasting controllers 
for each index are presented in section VII. 

Ill OPTIMAL INLET CONTROL 

The problem considered is not the complete inlet 
control problem with two degrees of control free- 
dom but rather is restricted to single-input, 
single -output control. This case is shown in 
Fig. 2 where the single control input is bypass 
door command and the one measurable output is 
throat exit static pressure (Pte) . This is a 
pressure measured Just downstream of the normal 
shock. The control is designed to take care of 
downstream disturbances only. Disturbances are 
considered to be random with a Gaussian probability 
distribution. 


Much of the work that has been done on inlet con- 
trol has been devoted towards modelling and simu- 
lation. Walittv 1 ) developed a Helmholtz resonator 
model for a mixed-compression inlet . Willoh' 2 / 
derived a linearized, small perturbation, high fre- 
quency model for an axi symmetric inlet. Mays (3) 
developed a nonlinear, large perturbation, low fre- 
quency model. All of the above approaches have 
been used in combination with test data to provide 
inlet control models. 

Controls development includes work done by Chun and 
BurrO 4 ) and Crosby, Neiner, and Cole. w) Designs 
for both of the aforementioned controls were ob- 
tained using frequency domain techniques. Barryw) 
conducted a design study based on an explicit des- 
cription of inlet disturbances. The disturbance he 
treated was atmospheric turbulence, described by - 
experimentally determined power spectral densities 
and probability distributions. Barry* s criterion 
for evaluating inlet controls was the expected fre- 
quency of unstart. This criterion, borrowed from 
the field of structural design, proved quite ap- 
plicable to inlet control and provided a meaningful 
criterion for Judging control effectiveness. 
Basically, it gives the designer a more rational 
basis (as opposed to frequency response methods) of 
Judging inlet control systems. 


In Fig. 2, the G*s represent transfer functions for 
the system dynamics, linearized about an operating 
point. GcoNTRGL the optimal controller, de- 
signed to control y s (shock position). Note that 
y s cannot be directly measured. Instead, the 
pressure P^ e is used. P te has been found to be 
related to shock position through shock position 
dynamics Gshqck. The additional transfer func- 
tions on Fig. 2 represent bypass door dynamics 
G3PP and noise coloring GCOLOR* Noise coloring 
accounts for the fact that the random disturbance 
at the compressor face is not in general white 
noise. Also shown in the figure is measurement 
noise on the signal Pt e * Ibis measurement noise 
is assumed white. 

A design constraint selected at the outset was that 
the control should be linear and time invariant. 
This is desirable first of all for simplicity. 
Secondly, linear time invariant control provides 
consistency to the analysis (Gaussian signals 
through linear systems remain Gaussian) . Also, and 
most important, linear optimal control theory re- 
sults can be used. 

The first of the design approaches taken is to 
minimize a quadratic performance index. This index 
can be written as 


In a previous work, (7) the authors extended Barry* s 
work to the design of optimal inlet controls which 
minimize the expected frequency of unstarts. This 
paper extends that work to a more realistic case 
in which actuator dynamics and nonwhite disturb- 
ances are included, as well as a comparison of al- 
ternate performance indices. Using a linearized 
operating point model, the present study compares 
and contrasts controls resulting from the minimiza- 
tion of two performance indices. The first index 
is the mean square value of shock position. The 
second is the expected frequency of inlet unstarts 
(as was used in ref. 6). 
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where the cr 2 terms are mean square values and k 
is a constant weighting parameter. Minimizing Jq 
minimizes the weighted sum of mean square shock 
position and bypass door command. The problem of 
minimizing Jq for a linear system is a classical 
linear quadratic problem. 


For the second approach a nonquadratic performance 
index was chosen to be minimized. This was given 
as 
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where 


Jnq. = * + < 2 l 
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Lambda is the expected frequency of unstarts, a is 
the shock position tolerance, and a 2, s are mean 
square values, This is a classic exceedance equa- 
tion and was used by Barry. ($) Due to the non- 
quadratic nature of J n q the standard linear quad- 
ratic results cannot be used directly. In this 
paper a desired linear time invariant control mini- 
mizing J n q was found using the quadratic equiva- 
lence principle. This technique was developed by 
Skelton^ 9 ) and applied to launch booster control 
system design. Use of the principle will be dis- 
cussed in more detail in section IV. Quadratic 
equivalence allows use of the same computer pro- 
grams which are used in solving the problem of 
minimizing Jq. The result obtained is the best 
linear control which causes the expected frequency 
of ms tarts to be a minimum for a given amount of 
control effort. 

TV OPTIMAL CONTROL THEORY 

The approaches to the inlet control design problem 
just described make use of the known results of the 
linear stochastic optimal control and estimation 
problem (Bryson and Ho, ref. 8) . The problem and 
solution can be summarized as follows; 

Given a time invariant plant described by 

x = Ax + Bu + Dw (3) 

with measurement vector 

z = Hx + v (4) 

and output vector 

y = Cx (5) 

x is the state vector, u the control vector and 
w and v are white Gaussian plant and measurement 
noise vectors. Plant and measurement noise cor- 
relation matrices are given by 

E(w(t)w T (t + t) ) = Q&(t) 

E{v(t)v^(t + t)} = R6(t) 

The problem is to minimize the quadratic perform- 
ance index 

J - E[1/2(x^Q c x + 2x?Nu + u T P e u)} (7) 

Equations (3) to (7) define the time invariant form 
of the linear stochastic optimal control and esti- 
mation problem. The solution is 

u := -Kc* (8) 

where the state estimate x is generated by a 




Kalman filter whose differential equation is 

k - Ax + Bu + Ke(z - H£) (9) 

Control gain matrix JQ, is 

Kc = Pc 1 (B T S + N T ) (10) 

where S Is the solution to the steady-state 
Biccati equation 

S(A - BP-1kT) + (A - BP“% T ) T S - S(BPq^B^)S 

+ CQc - NPJ 1 #) = 0 (11) 

Kalman filter gains are given by 

Ke = PH?B _1 (12) 

where P is the covariance of the estimation error 
and is the solution to 

AP + PA* - P(H T E _1 H)P + DQD T = 0 (13) 

If the error in the estimate is defined as 

e = £ - x (14) 

then the covariance matrices of x, % and e can 
be related by 

X = E(xx t ) = E{(£ - e)(£ - e) T ) = X + P (15) 

As shovm in Bryson, the covariance of % X, is 
given by the solution to the following Lyapunov 
equation 


(A - BKq)X + X(A - BKc)^ + P(H t R‘ 1 H)P = 0 (16) 

By adding Eqs. (13) and (16), a similar Lyapunov 
equation, 

(A - BK C )X + X(A - BKo) 1 + (BK C )P + P(BK 0 ) T 

+ DQD t = 0 (17) 

is obtained, which can be solved for the state 
covariance matrix X. Control gains Ke and 
Kalman filter gains Ke define the optimal con- 
troller. The state covariance matrix X contains 
mean square state information, which Is needed for 
overall system evaluation. 

Inlet Problem - Quadratic Performance Index 

Turning now to the particular case of interest, the 
inlet problem, the following scalar quantities are 
defined; u = u-^, w = w w , and y = y s . For the 
quadratic problem, it is desired that Jq(eq. (l)) 
be minimized. It is a simple matter to write Jq 
in form of the general quadratic index of Eq. (7). 
Since 

cr| s = E(yf) = E{x T (C T C)x) (18) 

and 
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= E{ku§), setting 

Q e = C T C 
N = 0 
P c = k 

allows the quadratic inlet problem to be solved 
using the solution given by Eqs. (8) to (17). 

Inlet Problem - Nonquadratic Performance Index 

Whereas the quadratic case of minimizing shock 
position can be solved in a straight forward 
manner, minimizing J n q (eq. (2)) requires some 
additional considerations. In this paper the quad- 
ratic equivalence principle will be used to convert 
the nonquadratic problem Into an equivalent quad- 
ratic one. For a discussion of quadratic equiva- 
lence, the reader is referred to Ref. 9. 

The necessary condition for existence of a sta- 
tionary point of J n q is that S(J n q) - 0 • Thus 
Eq. (2) becomes 

6(J n q) = 8(A) + kB^) (19) 

Expanding Eq. (19) in terms of cr 2 , oS , and cr 2 , 
one obtains Js ° 

0 = S(a| g ) + W x &(4 s ) + W 2 6(a| b ) (20) 

where parameters and Wg are given by 



The quadratic equivalence principle states that two 
performance Indices (one quadratic, one nonquad- 
ratic) are equivalent if their first variations are 
equal. This leads to the conclusion that the con- 
trol designed via quadratic equivalence is the best 
linear control that minimizes the nonquadratic per- 
formance index. 

In order to use the quadratic equivalence principle 
in the case under discussion, the variation of the 
quadratic performance index, J (eq. (t)), is set 
equal to zero, giving: 

8(E(x T Q c x]) + 2B(E(x T Nu) ) + b(E(u E 'p c u)) = 0 (22) 

To determine the conditions under which 5(J) = 
b(J n ^), it is necessary to expand Eq. (20) in 
terms of expected values. Since y s = CAx + CBu-^ + 
ClW-b, a|. g can be written as 



= E(y|y s } = E(CAx + CBu b + CDw w ) T (CAx 


+CBu t +CDw w )} (23) 

But C*D must equal zero else white noise w w 
would feed through to output y s such that o? 
would be infinite (making minimization of J n q y 
meaningless) . This is equivalent to saying that 
the transfer function relating w w to y g must 
have more poles than zeros. Expanding (23), using 
&y s from Eq. (18), then substituting for Oy s and 

crE in Eq. (20) and collecting terms, results in 

JS 

6(E(x T (A T C T CA + WjT^Cj.x}) + 8(E(2x T (A T C T CB)u t }) 


+ &(E{u^(B T G T CB + Wg)!^}) = 0 (24) 


Row compare Eq. (24), term by term, to Eq. (22). 
Solving the quadratic problem with 


Q 0 = A T C T CA + W X C T C 
IT = A T C T CB f 


(25) 


P c = b t c t cb + w 2 


gives the desired linear controller ( Kalman filter 
plus state estimate feedback) with one additional 
stipulation: At the optimum, Eqs. (21) must be 

satisfied. Since parameters and W2 are 

functions of mean square quantities which are in 
turn functions of the optimal control, an itera- . 
tive procedure is dictated. That is, select trial 
values of Wp and Wg and solve the quadratic 
problem using weighting matrices calculated using 
Eqs. ('25). Then using the quadratic solution (in 
particular, X, calculated from eq. (17)), new 
values of and W 2 can calculated. If 

these calculated values differ from the trial val- 
ues, a new trial pair is selected and the quad- 
ratic problem resolved. Particular iterative tech- 
niques are discussed in Ref. 9. For the results to 
be presented In this paper, however, simply re- 
running the quadratic problem for a set of (W^Wg) 
pairs proved adequate. 


One additional computation is required to calcu- 
late a £ and Oy from X, the state covariance 

matrix which is output from the quadratic problem. 
Since y s is a scalar it can be easily shown that 


a 2 = CXC T (26) 

y s 


and 

a? = E{y yT) = CAXA?C T - 2CA(X - P)kTbTCT 
y R s s c 


+ CBK c (X - P)k|b T C T (27) 

Also, the required mean square control effort can 
be expressed as 


<4 b = - p )^ (28) 
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This completes the formulation of the quadratic and 
nonquadratic inlet problems. The next section 
discusses the computer algorithms used. 

V COMPUTATION 

The first step in the computation procedure was to 
solve the estimation problem, yielding a set of 
Kalman filter gains Kg. This step is common to 
both the quadratic and nonquadratic problems, as 
both use the same Kalman filter to generate the re- 
quired state estimates. 

The nonquadratic problem was handled by solving the 
equivalent quadratic problem for a number of trial 
(Wi,W2') pairs. The solution to each quadratic 
problem K c was used with the Kg previously de- 
termined to compute the mean square values of the 
state vector. From this, Cy s , a|- s , and were 

calculated. For a range of k’ s, J n q, A, Wq, and 
Wg were then calculated. All these data were 
scanned off-line to determine optimum points. This 
was done by first finding points of minimum J n q 
for each value of k. As a cheek, an alternate 
method was used. This was to find points for con- 
stant k where the assumed Wq was sufficiently 
close to the calculated Wj and the assumed W3 
was close to the calculated ¥3 . The accuracy with 
which the Wq and W 2 equations (eqs. (21)) were 
satisfied was a function of the number of trial 
pairs (Wq,W 2 ). That is, finding the optimum more 
accurately required more pairs. The quadratic 
problem was solved like the nonquadratic problem 
except only one control problem solution was re- 
quired ( for a given value of k) . 

The two major subroutines used were one for the 
steady- state Riccati equation and one for the 
steady- state covariance matrix (Lyapunov) equa- 
tion. The steady- state Riccati equation subroutine 
used the negative exponential method. (10 ) The 
state covariance matrix equation (eq. (17)) was 
solved by transforming It into a set of n(n + l)/2 
linear equations. The solution method ehosen(H) 
was only one of a number of possible alternate 
methods for solving the Lyapunov equation. 

VI TRANSFER FUNCTIONS AND STATE VARIABLE MODEL 

Table I shows the various transfer functions re- 
quired to describe the inlet as shown in Fig. 2. 
^INLET and g SH 0CK transfer functions were 
obtained by curve fitting unpublished NA.SA data. 

The data were taken in wind tunnel tests on the 
NASA 40/60 axisymmetric inlet. Data for Crgpp 
appear in Ref. 12. Note that GjuleT an<3 - OSHOCK 
contain dead times. These dead times were ap- 
proximated by Fade approximations in order to ob- 
tain a finite state representation. The last item 
in the table, GCOLOR? represents the available 
knowledge concerning the compressor face disturb- 
ance. For example, if the disturbance were gener- 
ated at the upstream end, GgoLOR would serve to 
model the spectrum of atmospheric turbulence. 
However, for the downstream disturbance, the pole 
zero locations shown were chosen so that the color- 
ed disturbance is rich at low frequency (approxi- 
mating a step change in engine speed)., moderate at 
mid- frequencies (approximating combustion noise), 
and then cut off at high frequency. The disturb- 


ance spectrum of the output of this coloring cir- 
cuit is arbitrary but believed to be typical of 
that existing at an engine compressor face. 

Table II contains the necessary matrices which are 
used in the state variable description. These 
matrices were obtained using the transfer functions 
of table I. It can be seen in the A matrix that 
the first two states represent shock motion and the 
next four states represent inlet dynamics. The 
next two states are the noise coloring and the last 
two the bypass door dynamics. The value of measure- 
ment noise psd (r) was chosen from experience to 
approximate the level of noise on the throat exit 
static throat signal (Pte) found in wind tunnel 
tests. The psd of the white noise signal which 
drives the coloring circuit (q) is also arbitrary 
but was chosen so that the rms value of w<j is 
reasonable ( 0.168 kg/ sec) compared to a nominal in- 
let through flow of 16.8 kg/ sec. 

VII RESULTS 

As a basis for comparing both types of control, A -1 
(mean time between unstarts) was chosen as most 
representative of overall system performance. 

Since, in the performance index, the control vari- 
able was also weighted, results are shown with 
a wh (bypass door flow) as a parameter. In Fig 3, 

A-I is plotted as a function of 0^ with a 
(shock position tolerance) as the parameter. The 
solid curves are for the quadratic case, and the 
dotted curves are for the nonquadratic. Both de- 
signs use identical Kalman filters. The study used 
only one Kalman filter, since only one Q/R ratio 
was considered. 

The most significant point in Fig. 3 is that A -1 
is very sensitive to a. The value of A -1 in- 
creases with both a and (%,. For each case, 
each value of cw corresponds to a unique set of 
control gains ( j, Independent of a. It Is ap- 
parent in Fig. 3 that both cases are quite similar 
In their ability to prevent unstarts for a given 
value of control effort. The two curves for the 
nonquadratic case show that this control is only 
slightly better than quadratic control. The basic 
difference between the two cases is that in the 
nonquadratic case cr? is weighted in addition to 

a Js' 

For each value of a. A" 1 is bounded, for both 
large and small values of 0^, The system be- 
comes open loop as goes to zero. The open 

loop values of A"- 1 - could not be shown on Fig. 3 
because they are so small (for example, the open 
loop value of A~ 1 is 1.5x10“ hours for a = 

3.44 cm). As 0^ goes to infinity, for each 
value of a, A“^ doesn' t go to infinity (as it 
would were there no measurement noise), but instead 
reaches a limit. The limiting values of A" 1 
could not be conveniently computed, so are not 
shown on Fig. 3. The best that the control can do 
for large is to drive the estimate of x (not 

x. Itself) to zero. 

Of the many physical constraints present in the in- 
let control problem, the constraint on cr^ is 
most important. This is because the value of 
dictates the capacity of the bypass doors. 
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This, in turn, determines the size of the bypass 
door actuator. A practical question is, how large 
can <w be without incurring saturation on the 
bypass doors? 

Once the actuator is selected, the other con- 
straints must be considered. For the actuator used 
in these designs (described in ref. 12) there exist 
hard limits on bypass door velocity and acceler- 
ation. Thus, as an additional check, ms bypass 
door acceleration, position, and velocity were 
calculated. In determining how large various a 
values can be without incurring saturation, the 
following rule-of-thumb was used; make the three 
a value less than or equal to the limit value in 
question. Since the probability of exceeding the 
3cr limit is quite small, it is then unlikely that 
saturation will occur. For the quadratic case the 
use of the limits is illustrated as follows, using 
limit data from Ref. 12. 

Three critical values of <% are shown on Fig. 3. 
Point G is the value of rms bypass door flow at 
which the ms bypass, door position is equal to one 
third of its maximum (full open) value. If no 
other constraints existed, feasible optimal con- 
trollers would be all those producing ms bypass 
door flow rates less than or equal to 5.6 kg/ sec 
(point C) , However, when velocity and accelera- 
tion constraints on the bypass doors are consider- 
ed, the only feasible controllers are those which 
have rms bypass door flow rates less than or equal 
to 0.091 kg/sec (point A). Point A indicates the 
value of (and corresponding optimal controller) 

for which the rms bypass door acceleration is 
equal to one-third of its hard limit (maximum) 
value. Note that the velocity limit (point B) lies 
between points A and C, thus it is not the deter- 
mining factor in controller selection. A more 
general approach to the actuator saturation problem 
would have been to weigh, in the performance index, 
in addition to bypass door input, bypass door flow 
rate and rate of change of flow rate (both of which 
are state variables). This, however, was not done 
in this study. 

Controls are also compared on a closed loop fre- 
quency response basis. Figure 4 shows closed loop 
frequency responses for shock-position-to- 
compressor face flow disturbance. Note again the 
similarity between quadratic and nonquadratic de- 
signs. The fact that the colored noise w^ had an 
appreciable amplitude in the mid frequency range 
caused the closed loop to stay under the open loop 
until about 70 Hz. Zero frequency disturbances are 
attenuated by about two orders of magnitude. Re- 
ferring to Fig. 3, for the levels of control effort 
selected, it can be seen that small changes in the 
frequency response magnitude (fig. 4) correspond to 
large changes in A' 1 . This again demonstrates the 
difficulty in using A alone as the performance 
criterion - it is very sensitive to small changes 
in controller parameters. 

In addition to closed loop analyses, the Kalman 
filter -plus -control gain combinations were examined 
as conventional controllers. Figure 5 shows the 
open loop frequency response of bypass door com- 
mand voltage to a change in throat exit static 
pressure (transfer function ^CONTROL in 2 ) • 


This corresponds to the feedback compensator of 
classical control design. Again, quadratic de- 
signs are seen to be very similar to nonquadratic 
designs (for the same value of rms bypass flow). It 
should be noted here that the inlet is uncontrolla- 
ble, but that the uncontrollability creates no 
problem with any of the computer algorithms. In 
particular, it can be seen in table I that the 500 
rad/ sec pole In GsHOCK is ‘uncontrollable, as are 
the two noise coloring poles. Since all are stable, 
the uneontr pliability creates no problems. 

The controller represented in Fig. 5 is tenth 
order. Since a tenth order control is possibly 
more complex than desired, simplification techni- 
ques were investigated. One method of simplifica- 
tion would be to curve fit the data on Fig. 5 with 
a model of desired order. A better approach to 
simplification would be to apply a general simpli- 
fication technique (for example, as shown in 
ref. 13) to the controller. For this problem, it 
turned out that basic simplifications could be 
made Just by examining the controller root locus. 

Two pole zero cancellations and one near cancella- 
tion can be seen on Fig. 6. Figure 6 was plotted 
as a root locus of the controllers for the quad- 
ratic case. It was made by taking the eigenvalues 
of (A - BKq - KgH) and zeros of Kg (si - (A - BK,, - 
KgH) ) . The parameter on the loci is cr™ . Of 

the ten poles and nine zeros shown, three cancel- 
lations occur. In particular, the two thousand 
rad/ sec noise coloring pole is almost cancelled 
in addition to the 800 rad/sec and 500 rad/sec 
shock position poles which are exactly cancelled. 
Thus, the controller can be implemented as a 6th/ 
7th order system. It's possible, of course, that 
further simplifications could be made and results 
of the simplifications evaluated by recomputing Jg. 

VIII CONCLUSIONS 

It was found that through the use of the quadratic 
equivalence principle, a control was designed which 
minimized a physically meaningful performance index. 
The control system resulting is the best linear 
control which minimizes the nonquadratic index. 

The particular nonquadratic index used was quite 
sensitive to shock position tolerance. It was 
found that using a quadratic index (minimizing 
shock position deviation) gave results very simi- 
lar to the nonquadratic results. Since the quad- 
ratic approach requires only one computation of the 
Riccati equation for each controller, while the 
nonquadratic requires several, the simpler quad- 
ratic approach may still be more attractive. While 
the example in this study was a single -input, 
single -output case, the methods could be extended 
to the more realistic inlet problem, including 
atmospheric turbulence disturbance, centerbody 
actuation for throat Mach number control, and com- 
pressor face or other additional pressure measure- 
ments. 
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♦Blanks are all zeros. 

q - psdCwb) = 0.188 (kg/sec) 2 /rad/sec 

r = psd(vp^) = 4.754X10 -4 ( N/ cm 2 ) / r ad/se c 


TABLE XI. - HjANT MATRICES AND NOISE PSD*S 
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Oblique shock 



Figure 1. - Schematic of su 



OPTIMAL CONTROLLER - KALMAN 
FILTER PLUS STATE ESTIMATE 
FEEDBACK GAINS 


Figure 2. - Inlet control system block diagram. 









-6512 


QUADRATIC CASE 

— — NONQUADRATIC CASE 



Figure 4 - Closed loop frequency responses, quadratic and nonquadratic 
cases, shock position to compressor face flow disturbance. 



FREQUENCY, Hz 

Figure 5. - Controller frequency responses, quadratic, and non-" 
quadratic cases, bypass door command voltage to throat exit 
static pressure. 




Figure 6, - Controller root locus for quadratic case with o w asthe 
parameter. b 
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